Nematic phase of the two-dimensional electron gas in a magnetic field 
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The two dimensional electron gas (2DEG) in moderate magnetic fields in ultra-clean AlAs-GaAs 
heterojunctions exhibits transport anomalies suggestive of a compressible, anisotropic metallic state. 
Using scaling arguments and Monte Carlo simulations, we develop an order parameter theory of 
an electron nematic phase. The observed temperature dependence of the resistivity anisotropy 
behaves like the orientational order parameter if the transition to the nematic state occurs at a finite 
temperature, T c ~ 65mK, and is slightly rounded by a small background microscopic anisotropy. 
We propose a light scattering experiment to measure the critical susceptibility. 
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Recently^ two of us introduced the concept of liq- 
uid crystal phases of the two-dimensional electron gas 
(2DEG) in large magnetic fields, as an extention-,of ear- 
lier work on high temperature superconductors&tl Elec- 
tronic liquid crystal phases are quantum mechanical-ana- 
logues of classical liquid crystals, and are predictedcm to 
be a generic feature of strongly correlated fermionic sys- 
tems. In the case of the 2DEG the competing effects 
of repulsive (Coulomb) interactions and the quenching 
of the kinetic energy of electrons in Landau levels lead 
naturally to the existence of such phases. Pursuing this 
analogy, the phases of a 2DEG in order of increasing 
symmetry breaking, were characterized as: (a) isotropic 
liquids, (b) nematic liquids, (c) smectic liquid crystals, 
and (d) insulating crystals. The fluid character of states 
(a) and (b) is obvious, as they are translationally invari- 
ant. The smectic breaks transational symmetry in only 
one direction, and so is also a fluid. Insulating crystals 
break the translational symmetry down to a discrete sub- 
group, such that there are an integer number of electrons 
per unit cell. 

Among the isotropic liquids are the various quan- 
tum Hall fluids, while insulating crystals are simply 
the Wigner crystal and its generalizations. A smec- 
tic (or stripe) phase has been found in Hartree-Fock 
calculations, presumably accurate in high Landau lev- 
elsaill, which provides a qualitative picture of a smectic 
state. Moreover, a recent exact diagonalization study 
of a system with 12 electrons in the N > 2 Landau 
levels found -results consistent with a smectic (stripe) 
ground stateQ (up to yet poorly understood finite size ef- 
fects). Wigner crystals and "bubble phases" (crystalline 
states with several particles per unit cell) have also been 
found as Hartree-Fock variational states. In addition, a 
predicted charge-density-waifeJpstability of the smectic 
phase at Hartree-Fock levelEJoErEiJ leads to an insulating 
"stripe-crystal" phase with a rectangular unit cell. How- 
ever, a microscopic theory of the nematic phase does not 
presently exist. 

The low energy physics of a quantum smectic can be 
understood in terms.xi£-.a theory of quantum fluctua- 
tions of the smecticD'QElll3, including the phase transition 



to the stripe crystal phase. The long-distance behavior 
of a quantum nematic phase is completely determined 
by symmetry and by the associated Goldstone modes. 
Both smectic and nematic states break rotational sym- 
metry, and as such transport properties of both states is 
anisotropic. In Ref. pj we argued that while the smectic 
and nematic states are both natural candidates to- eps-, 
plain the anisotropy observed in recent experimentaiirtHI 
on 2DEG in high mobility AlAs—GaAs heterostructures, 
at least at finite temperature, there are strong reasons to 
favor the nematic. The smectic has an infinite conduc- 
tivity, at least in one direction, whereas the measured 
(anisotropic) conductivity has a finite T — > limit. In 
addition, the data shows a pronounced temperature de- 
pendence of the resistivity, consistent with the existence 
of a finite temperature phase transition; since the energy 
of a dislocation is still finite even for a Coulomb inter- 
actioH-thc smectic always melts at any non-zero temper- 
atureEJ. In the present paper we explore the universal 
properties of the 2D electron nematic. 

The experiments of Refs. O and 03 have revealed 
the existence of regimes of magnetic fields in which the 
2DEG exhibits characteristics of a compressible fluid 
with an unexpectedly large and temperature-dependent 
anisotropy in its transport properties. This occurs when 
the Landau level index N lies in the range 2 < N < 6. 
The large anisotropy is seen only in high mobility sam- 
ples. In the same samples, similar behavior has also been 
seen in the first N—r= 1 Landau level when the mag- 
netic field is tiltedli-l while in the lowest N = Lan- 
dau level, a large number of fractional quantum Hall 
(FQH) states are observed, including a subtle FQH state 
at v = 5/2. To summarize the experimental facts: (i) 
Measurements on square samples show that, as the tem- 
perature is lowered below 100 mK, the longitudinal resis- 
tance R xx grows very rapidly while R yy becomes smaller; 
as T — > 0, their ratio approaches a constantli3 in the range 
100 < Rxx/ Ryy < 3, 500 where aj-apd y are orthogonal 
lattice directions, (ii) In Hall baraliO, R yy is essentially 
temperature independent while R xx increases by a factor 
of 5— 10 as the temperature is lowered below IQQmK . (iii) 
The compressible (dissipative) regime occupies a finite 



range of magnetic fields AB, centered around the middle 
of the partially filled Landau level, and unlike the conven- 
tional transition between plateaus, AB does not shrink 
as the temperature is lowered, (iv) In the same range 
of magnetic fields the Hall resistance varies continuously 
with the magnetic field, (v) At the peak, the conductiv- 
ity a yv is typically of the order of e 2 /h (see below), (vi) 
Applying an in-plane magnetic field tunes the anisotropy 
through zero, and reverses the roles of the x and y direc- 
tions, (viii) In the region of the resistivity peak, the I— V 
curves are highly non-linea* but do not show any thresh- 
old (depinning) behavior Ota. (ix) There are reentrant 
integer quantum Hall plateaus symmetrically located for 
magnetic fields outside the compressible regime, (x) The 
anisotropy has not been reported in lower mobility sam- 
ples which show instead the (expected) phase transition 
between quantum Hall states. 

A natural interpretation of the experiments is that, 
in regimes in which interactions dominate over the ef- 
fects of disorder, instead of the expected transition be- 
tween plateaus in the middle of the partially filled Landau 
level, the 2DEG forms a compressible anisotropic fluid. 
Because a continuous (rotational) symmetry cannot be 
spontaneously bppken in D=2, for such an anisotropy 
to be observable^ the sample must have a small back- 
ground microscopic anisotropy whose effect is greatly 
amplified at low temperatures by the collective prop- 
erties of the state. Specifically, we will show that the 
experimentally observed temperature dependence of the 
anisotropy can be understood as evidence for a finite tem- 
perature Kosterlitz-Thouless (KT) transitipixfrom a two- 
dimensional nematic to an isotropic fluidtjO, which is 
rounded by a symmetry breaking field representing the 
effects of the background anisotropy. We present a&,anal- 
ysis (Fig. 1) of the experimental data of Lilly et aO, and 
a fit with the results of a Monte Carlo simulation of a 
model of a classical nematic in a symmetry breaking field 
(Fig. 2). The results strongly support our earlier claimEJ 
that the anisotropic transport occurs where the 2DEG is 
in a nematic phase (at least at finite temperature). They 
also give some indirect support to the further conjecture 
made in Ref. [jL that there is a direct transition as a func- 
tion of B from the nematic state to an insulating stripe 
crystal phase, which was identified with the innermost of 
the reentrant quantum Hall liquids. (See, also Ref. [18].) 

In 2D, the nematic has only quasi- long-range order; its 
finite temperature transition to a disordered liquid can 
be described by the two-dimensioaal classical AF-model 
with a director order parametera. Such a description 
should fail at (very) low temperatures where quantum 
fluctuations (and/or quenched disorder) become impor- 
tant. Since the order parameter of the nematic state is 
a director field, m(r), it is periodic under rotations by 
7T, and has the form m x (r) + im y (r) = exp(2i9p). The 
classical Hamiltonian of this system is thus 



H = -J ^2 cos(2A M 6» r -) + h^2 cos(26> ? ) 



(1) 



where, for simplicity we have used a square lattice of unit 
spacing whose sites are labelled by the lattice vectors r. 
In Eq. [fl we have use the notation A^9p — Op + g — Op, 
where e^ is a unit vector along the direction /i = x,y, 
and J is the stiffness, the energy required to rotate two 
nearby regions by a small angle. The quantity h breaks 
rotational symmetry explicitly. It represents the effects 
of a background symmetry breaking field, such as the 
anisotropy and/or the effects of a parallel magnetic field. 

Because two is the lower critical dimension for contin- 
uous symmetry breaking, for < Tp< T c the system is 
controlled by a line of critical pointai3. In this range of 
T, where classical Goldstone excitations of the nematic 
("spin waves") dominate, the correlation function G{r) = 
(exp(2i[#(r) — 0(0)]) of the order parameter has power 
law behavior, G{f) ~ l/\r\^ T \ with n = 2T/ttk(T) and 
a divergent susceptibility. Here n — ► 1/4 as T — > T c and 
k(T) is the helicity modulus, which approaches k(0) = 4 J 
as T — > and n(T c ) — 4T c /tt. In the presence of a 
symmetry breaking field h, the order parameter behaves 
like to = (exp(2i#)) ~ l^l 1 ^ 15 , where 5 — 4/ry — I and 
S(T C ) — 15. For T > T c , the correlation length is finite, 
and diverges at T c like £ ~ exp(A/^/T — T c ), where A 
is a (non-universal) constant. At finite h the correlation 
length is always finite, and the singularities of the KT 
transition get rounded. In this regime, even a very small 
symmetry breaking field induces a very large expectation 
value of the order parameter. For T > T c , to ~ x(T)/i, 
where x(T) ~ £ 7 / 4 is the susceptibility. 

How is this thermodynamic transition, which describes 
the breaking of rotational invariance, related to the trans- 
port properties? On general grounds one expects that 
near a phase transition, quantities which transform in 
the same way under the symmetry should be related, 
even if one is a transport coefficient .and the other a ther- 
modynamic property. In particulaiu the combination of 
resistivities £ = (p xx — p yy )/(p xx + Pyy) transforms like 
the order parameter m = (exp(i26)) . It should be related 
to the order parameter through an odd analytic function 
C = /(to). Therefore, near T c , if the symmetry breaking 
is small, the linear approximation /(to) oc to + 0(m 3 ) 
should be reasonably accurateB3. 

We can determine if the 2DEG is in a nematic phase 
by analyzing the temperature dependence of the resis- 
tivity in terms of the temperature dependence of the or- 
der parameter of the nematic in the presence of a sym- 
metry breaking field. What is needed is the function 
to = $(T, h), the equation of state, which we computed 
by a Monte Carlo simulation of the classical XY model of 
Eq. II]. Notice that we relate to to a local (intensive) prop- 
erty such as the resistivity instead of to the resistance, 
which is extensive and sensitive to significant finite size 
effects. However, the experimental data gives the resis- 
tances R xx and R yy as functions of temperature, not the 
resistivities. Thus, in order to fit the data, we extracted 
the resistivities from the measured resistances, using a 
method discussed below, with the result shown in Fig. [P. 
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FIG. 1. Resistivities p xx and predetermined from the resis- 
tance data of Lilly et al 111] fl6| at v = 9/2, after deconvolv- 
ing the effects of the geometry; p yy is essentially constant for 
the entire range of temperatures, as in Hall bars. Inset: the 
function x(T). 

One result of this analysis is that the resistivities for 
the square sample behave exactly in the same way as the 
resistances of the Hall bars. Given the low T values of 
p xx and pyy in Fig. 1, and the (large) measured value 
of the Hall conductance, one finds that the peak value of 
the conductivity is a yy — 1.12e 2 /h and a xx = 0.11e 2 /h. 
Notice that p xx saturates rather sharply below bbmK 
and that both p xx and p yy approach non-zero (and dif- 
ferent) values as T — > 0. Thus, the 2DEG remains in 
an anisotropic compressible (metallic) state, down to the 
lowest accessible temperatures. 
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FIG. 2. Fit of the Monte Carlo data for a 100 x 100 lattice, 
to the data of Lilly et. al @ @. The best fit is found for 
J = 73mK and h = 0.05J = 3.5mK and T c = 65mK . 

Having determined the temperature dependence of 
the resistivities we can now see if it is consistent with 
a (rounded) phase transition from a high temperature 
isotropic fluid phase to a low temperature nematic phase. 
We have done this by means of a Monte Carlo Metropolis 



simulation of the 2D XY model of Eq. [I] on square lat- 
tices of sizes 40 x 40 through 120 x 120, for the range 
of symmetry breaking fields 0.01J < h < 0.5 J, and 
for a wide range of temperatures (see below). In Fig. 
@ we show our Monte Carlo data for the order param- 
eter as a function of temperature for h = 0.05J. For 
this range of symmetry breaking fields, we find that for 
L = 100 the finite size effects on the order parameter are 
very small. We have fitted the data by assuming that 
C = {Pxx - PywhiiPxx + Pyy) is actually equal to the order 
parameter mrJ Having done so, we fitted the data by 
finding the best value of J that fits the data for a given 
value of h, and then changed h to get the best fit. 

The classical nematic does indeed explain the temper- 
ature dependence; the data is consistent with a thermo- 
dynamic KosterlitZj-Thouless transition at T c (h = 0) = 
0.88J ~ 65mK Ej, slightly rounded by a background 
anisotropy field of magnitude h ~ 0.05 J ~ 3.5mK, which 
is a very small energy scale. Notice that both the stiffness 
J ~ 73mK and h are much smaller than the Coulomb 
energy, although they are comparable with the gap in 
the v = 5/2 (presumably paired) state, which hints a 
possible common origin. Below 55771^ the fit is not as 
good. In this temperature range the XY order parame- 
ter is big (larger than 1/2) so there is no reason to expect 
£ ~ m. However, £ strikingly saturates (unlike m which 
shows the characteristic linear temperature behavior of 
classical spin wave theory) so the discrepancy may be 
indicating that quantum mechanical effects (or disorder) 
are important at low T. 
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FIG. 3. Longitudinal and transverse susceptibilties of a 
classical nematic on a 100 x 100 lattice, for h = 0.05J. 

Our analysis of the experiments strongly indicates 
that the 2DEG in large magnetic fields in clean sam- 
ples has regimes where it behaves like a nematic fluid, 
an anisotropic metal. Such a metallic state should have 
a strong signature in polarized light scattering experi- 
ments. In particular, a nematic has long range fluctu- 
ations in the orientational order, which will cause the 



polarization tensor correlation function (and the corre- 
sponding longitudinal and transverse susceptibilities xl 
and xt) to have a singularity at T c (cut off by the 
anisotropy). This effect is similar to critical opalescence 
but for orientational order instead of density fluctua- 
tions. For non-zero and small background anisotropy 
h, for T < T c , xl is XL ~ h~ a , with a = 1 — 1/5, 
where a = 14/15 at T c . For T > T c , xl can be writ- 
ten in a scaling form as x{h,T) ~ £ 7/4 $o(^£ 15 ^ 8 ) where 
$ (0) = 1 and $ (x) ~ a;^ 14 / 15 as x — > oo; as dis- 
cussed above, £(T) ~ exp(- J 4/v / t), where £ = T/T c - 1. 
Thus, at fixed but small h, as the temperature is low- 
ered, xl increases very rapidly to a maximum above 
T c , with a crossover to a critical behavior ~ h~ a ( T \ 
where a(T) = 2(7 + £)/(15 + i) for \h\ < T, and 
Xl~ (7rT/2/i)(l/(47r 2 J + /i)), forT< |/i|. On the other 
hand, by Goldstone's Theorem, xt = m/h, where m is 
the order parameter (Fig. 0). xl and xt are shown in 
Fig. g. Although the specific heat of the 2DEG is very 
hard to measure, it should have a broad bump at a tem- 
perature T* > T c set by the core energy of the disclina- 
tions of the nematic phase, while at T c there is a very 
weak essential singularitylill (rounded by the anisotropy) . 
Finally, we used the following procedure to extract the 
resistivities frqm.-thc measured resistances. It was ob- 
served recentlyOO that much of the discrepancy between 
the data from Hall bars and square samples can be un- 
derstood to be a consequence of a large distortion of the 
current distribution in square samples. If one assumes 
that the 2DEG_has an anisotropic but homogeneous re- 
sistivity tensorEj, one can calculate the distribution of 
currents using the method of conformal mappingscj. One 
finds that, for an anisotropic homogeneous sample, with 
its principal axes aligned with the edges, 



Rxx/Ryy = g(x)/g(l/x) 



(2) 



where x = \L y /L x \^J p yy / p xx = x(T) measures the as- 
pect ratio L y /L x and the ratio of resistivities, and g(x) 
is given by 



g(x) = ln( 



9 3 (inx/2) + Vk0 2 (iTrx/2) , 



9 3 (iirx/2) - Vk6 2 (iTrx/2) 
#2(2) and #3(2) are theta-functions with modulus k 



1 , „2n 



(3) 



(4) 



where q — exp(— 2irx) is the period c3. Given R xx /R yy at 
different temperatures, and using Eq. (0), we calculated 
the function x(T) (shown in the inset of Fig. |l|). At high 
T, x(T) approaches a value somewhat larger than 1, but 
it is smaller than 1 at lower temperatures, and both the 
resistances and the. resistivities show a crossing at some 
high temperatureo. This effect indicates that the sam- 
ple is not homogeneous at large scales. A macroscopic 
inhomogeneity is equivalent to an effective aspect ratio 



L y /L x 7^ 1, and by choosing a value for L y /L x = 1.12 
we can make the ratio p xx / ' p yy approach unity at high 
temperature. 
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